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Abstract 

Unambiguously distinguishing between nonorthogonal but linearly independent quan- 
tum states is a challenging problem in quantum information processing. In this work, an 
exact analytic solution to an optimum measurement problem involving an arbitrary num- 
ber of pure linearly independent quantum states is presented. To this end, the relevant 
semi-definite programming task is reduced to a linear programming one with a feasible 
region of polygon type which can be solved via simplex method. The strength of the 
method is illustrated through some explicit examples. Also using the close connection 
between the Lewenstein-Sanpera decomposition(LSD) and semi-definite programming 
approach, the optimal positive operator valued measure for some of the well-known ex- 
amples is obtain via Lewenstein-Sanpera decomposition method. 

Keywords: Optimal Unambiguous State Discrimination, Linear Programming, Lewenstein- 
Sanpera decomposition. 

PACs Index: 03.67.Hk, 03.65.Ta, 42.50.-p 
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In quantum information and quantum computing, the carrier of information is some quantum 
system and information is encoded in its state. A quantum state describes what we know 
about a quantum system. Given a single copy of a quantum system which can be prepared in 
several known quantum states, our aim is to determine in which state the system is. This can 
be well-understood in a communication context where only a single copy of the system is given 
and only a single shot-measurement is performed. This is in contrast with usual experiments 
in physics where many copies of a system are measured to get the probability distribution 
of the system. In quantum state discrimination, no statistics is built since only a single-shot 
measurement is performed on a single copy of the system. Actually there are fundamental 
limitations to the precision with which the state of the system can be determined with a single 
measurement. A fundamental aspect of quantum information theory is that non-orthogonal 
quantum states cannot be perfectly distinguished. Therefore, a central problem in quantum 
mechanics is to design measurements optimized to distinguish between a collection of non- 
orthogonal quantum states. The topic of quantum state discrimination was firmly established 
in 1970s by pioneering work of Helstrom [1], who considered a minimum error discrimination of 
two known quantum states. In this case, the state identification is probabilistic. Another pos- 
sible discrimination strategy is the so-called unambiguous state discrimination (USD) where 
the states are successfully identified with non-unit probability, but without error. USD was 
originally formulated and analyzed by Ivanovic, Dieks and Peres [2, 3, 4] in 1987. The solu- 
tion for unambiguous discrimination of two known pure states appearing with arbitrary prior 
probabilities was obtained by Jaeger and Shimony[5]. Although the two-state problem is well 
developed, the problem of unambiguous discrimination between multiple quantum states has 
received considerably less attention. The problem of discrimination among three nonorthogo- 
nal states was first considered by Peres and Terno [4] . They developed a geometric approach 
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and applied it numerically on several examples. A different method was considered by Duan 
and Guo [6] and Y. Sun and et al [7]. Chefles [8] showed that a necessary and sufficient con- 
dition for the existence of unambiguous measurements for distinguishing between quantum 
states is that the states are linearly independent. He also proposed a simple suboptimal mea- 
surement for unambiguous discrimination for which the probability of an inconclusive result is 
the same regardless of the state of the system. Equivalcntly, the measurement yields an equal 
probability of correct detection of each one of the ensemble states. 

Over the past years, semidefinite programming (SDP) has been recognized as a valuable 
numerical tool for control system analysis and design. In SDP, one minimizes a linear function 
subject to the constraint that an affine combination of symmetric matrices is positive semidef- 
inite. SDP has been studied (under various names) as far back as the 1940s. Subsequent 
research in semidefinite programming during the 1990s was driven by applications in combi- 
natorial optimization [9], communications and signal processing [10, 11, 12], and other areas 
of engineering [13] . Although semidefinite programming is designed to be applied in numerical 
methods, it can be used for analytic computations, too. In the context of quantum computa- 
tion and quantum information, Barnum, Saks and Szegedy have reformulated quantum query 
complexity and in terms of a semidefinite program [14], while M. A. Jafarizadeh, M. Mirzaee 
and M. Rezaee have revealed the close connection between Lewenstein-Sanpera decomposi- 
tion, robustness of entanglement, finite quantum tomography and semi-definite programming 
algorithm. [15, 16, 17, 18]. 

The problem of finding the optimal measurement to distinguish between a set of quantum 
states was first formulated as a semidefinite program in 1972 by Holevo, who gave optimality 
conditions equivalent to the complementary slackness conditions [1] . Recently, Eldar, Megret- 
ski and Verghese showed that the optimal measurements can be found efficiently by solving 
the dual followed by the use of hnear programming [19] . Also in [20] , SDP has been used to 
show that the standard algorithm implements the optimal set of measurements. All of the 
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above mentioned applications indicate that the method of SDP is very useful. The reason 
why the area has shown relatively slow progress until recently within the rapidly evolving 
field of quantum information is that it poses quite formidable mathematical challenges. Ex- 
cept for a handful of very special cases, no general exact solution has been available involving 
more than two arbitrary states and mostly numerical algorithms are proposed for finding op- 
timal measurements for quantum-state discrimination, where the theory of the semi-definite 
programming provides a simple check of the optimality of the numerically obtained results. 

In this study, we obtain the feasible region in terms of the inner product of the states 
which enables us to solve the problem analytically without using dual states. Exact analytical 
solution for optimal unambiguous discrimination of N linearly independent pure states is 
calculated and a geometrical interpretation for minimum inconclusive result for unambiguous 
discrimination of two pure states is presented. For more than three states, the analytical 
calculations is very complicated to write down and therefore we will consider the spacial cases 
such as geometrical uniform states and Welch bound equality (WBE) sequences. To solve the 
problem in general form, following prescription of Refs [21, 22, 23, 24, 25], we have reduced it 
to LP one, where the computation can be done at a very fast pace and with high precision. 
Moreover, we obtain the feasible region in terms of the inner product of the dual states and 
show that LSD is equivalent to optimal unambiguous state discrimination, and thus one can 
use LSD to solve the problem of optimal unambiguous state discrimination. This method is 
illustrated for two and three linearly independent states explicitly. 

The organization of the paper is as follows: First the definition of the unambiguous quan- 
tum state discrimination are given. After that in section 3 unambiguous discrimination of 
quantum states by introducing feasible region and using linear programming are discussed. 
Then Lewnstein-sanpera decomposition is studied as an optimal unambiguous discrimination 
of quantum states. Finally, discrimination of non-orthogonal quantum states using approxi- 
mated linear programming are discussed. The paper is ended with a brief conclusion and an 
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appendix. 

2 Unambiguous quantum state discrimination 

In quantum theory, measurements are represented by positive operator valued measures (POVMs). 
A measurement is described by a collection {Mk} of measurement operators. These operators 
are acting on the state space of the system being measured. The index k refers to the mea- 
surement outcomes that may occur in the experiment. In quantum information theory the 
measurement operators {M^.} are often called Kraus operators [26]. If we define the operator 

n, = MlMk, (2-1) 

the probability of obtaining the outcome k for a given state pi is given by p{k\i) — Tr{IlkPi)- 
Thus, the set of operators is sufficient to determine the measurement statistics. 

Definition (POVM). A set of operators {11^} is named a positive operator valued measure 
(POVM) if and only if the following two conditions are met: 

1. Each operator 11^ is positive (V' | 11^. | ■0) > 0, V | ■0)- 

2. The completeness relation is satisfied, i.e., 

Enik = l- (2-2) 

k 

The elements of {Ilfc} are called effects or POVM elements. On its own, a given POVM {Ilfc} 
is enough to give complete knowledge about the probabilities of all possible outcomes; mea- 
surement statistics is the only item of interest. Consider a set of known states pi,i — 1, N, 
with their prior probabilities rji. We are looking for a measurement that either identifies a state 
uniquely (conclusive result) or fails to identify it (inconclusive result). The goal is to minimize 
the probability of inconclusive result. The measurements involved are typically generalized 
measurements. A measurement described by a POVM {11^}^^ is called unambiguous state 
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discrimination measurement (USDM) on the set of states {pijfLi if and only if the following 
conditions are satisfied: 

1. The POVM contains the elements {Ilk}k=o where N is the number of different signals 
in the set of states .The element Ho is related to an inconclusive result, while the other 
elements correspond to an identification of one of the states pi, i — 1, N. 

2. No states are wrongly identified, that is, Tr{pin.k) —0 \/i i,k — 1, N. 

Each USD measurement gives rise to a failure probability, that is, the rate of inconclusive 
result. This can be calculated as 

g = ^r/,Tr(Tr(p,no)). (2-3) 

i 

And the success probability can be calculated as 

P = 1 - g = Y.ri{rr{Tr{p,Ili)). (2-4) 

i 

A measurement described by a POVM {H^^*} is called an optimal unambiguous state dis- 
crimination measurement (OptUSDM) on a set of states {pi\ with the corresponding prior 
probabilities {r]i\ if and only if the following conditions are satisfied: 

1. The POVM {n^P*} is a USD measurement on {pi} 

2. The probability of inconclusive result is minimal, that is, Q{{Xl°^^}) — min(5({nfe}) 
where the minimum is taken over all USDM. 

Unambiguous state discrimination is an error-free discrimination. This implies a strong con- 
straint on the measurement. Suppose that a quantum system is prepared in a pure quantum 
state drawn from a collection of given states < i < N} in d-dimensional complex 

Hilbert space H with d > N. These states span a subspace U oi H. In order to detect the 
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state of the system, a measurement is constructed comprising + 1 measurement operators 
{rij, < i < N}. Given that the state of the system is the probabihty of obtaining out- 
come k is {ipilUklipi) ■ Therefore, in order to ensure that each state is either correctly detected 
or an inconclusive result is obtained, we must have 

(iPilUklA) = PiSik, l<t,k<N, (2-5) 

for some < < 1. Since Uq — — Y,^i Hj, we have {ipi\Ilo\jJji) — 1 — pi. So a system 
with given state the state of the system is correctly detected with probability pi and 
an inconclusive result is obtained with probability 1 — Pi. It was shown in [8] that (2-5) is 
satisfied if and only if the vectors are linearly independent, or equivalently, dimW = dimH. 
Therefore, we will take this assumption throughout the paper. In this case, we may choose 
[27] 

Ui^Pi\i,i){il;i\, l<i<N, (2-6) 

where, the vectors e U are the reciprocal states associated with the states i.e., there 
are unique vectors in U such that 

{AM^Sik, l<i,k<N. (2-7) 

With $ and l» we denote the matrices such that their columns are j'^j) and ji/'j), respectively. 
Then, one can show that $ is 

|) = $($$*)-!. (2-8) 

Since the vectors i = 1, N are linearly independent, $$* is always invertible. Alterna- 
tively, 

I, ^ ($$*)t$^ (2-9) 

so that 

l^i) = ($r)*|^,) (2-10) 
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where, (.)■'• denotes the Moore- Penrose pseudo-inverse [28]. The inverse is taken on the subspace 
spanned by the columns of the matrix. If the state is prepared with prior probabihty rji, 
then the total probabihty of correctly detecting the state is 

N N 

P = Y.Vrmiii\A) =T.V^P^ (2-11) 
i=l 1=1 

and the probability of the inconclusive result is given by 

TV N 

Q=l-P = J2v^{^^\^0\i^^) = 1 ' T.V^P^■ (2-12) 
1=1 i=l 

In general, an optimal measurement for a given strategy depends on the quantum states and 
the prior probabilities of their appearance. In the unambiguous discrimination for a given 
strategy and a given ensemble of states, the goal is to find a measurement which minimizes the 
inconclusive result. In fact, it is known that USD (of both pure and mixed states) is a convex 
optimization problem . Mathematically, this means that the quantity which is to be optimized 
as well as the constraints on the unknowns, are convex functions. Practically, this implies that 
the optimal solution can be computed in an extremely efficient way. This is therefore a very 
useful tool. Nevertheless our aim is to understand the structure of USD in order to relate 
it with neat and relevant quantities and to find feasible region for numerical and analytical 
solutions. So, by using SDP we determine feasible region via reciprocal states and reduce the 
theory of the SDP to a linear programming one with a feasible region of polygon type which 
can be solved via simplex method exactly or approximately. 



3 Unambiguous discrimination of quantum states using 
linear programming 

The method presented in this section, seems to be a powerful method which enables us to 
analytically discriminate N linearly independent pure quantum states. We naturally want to 
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the probabilities pi be as large as possible in order to increase the detection probabilities, but 
their values are bounded by the demand of positivity of Ho. Let 



N 



i=l 



with {ipi I ip-^) = 0, and normalization condition is defined by 



TV 



I ^) = ^ a>,-(V^, I + (^^ I = 1. 



(3-13) 



The relation (3-13) can be written as 



AT 



J2 (AajGij < 1, 



(3-14) 



where, Gij — {ipi \ ipj) are matrix elements of the Gram matrix. Now we define the following 
vector representation 



ai 



X 



(3-15) 



then, the Eq. (3-14) can be rewritten as the following constraint 

X^GX < 1. 



(3-16) 



Positivity of Tr (IIo) gives 

i 

This last condition is a decisive one that actually determines the domain of acceptable values 
of Pi. This result leads us to the optimization problem defined as 

maximize J2iViPi ^ 1 

j:ia^Pi<l (3-17) 
X^GX < 1. 

If we write pi — A'^j with < < 1, then we will have 

Y: a^Pi = ^Ui < 1 ^ A' < -1^. (3-18) 



S.t 
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Then, we must compute the maximum value of J^i (ij^i such that A' possesses its lowest possible 
value. That is, 

maximize J2i c^^i 
s.t X^GX < 1. 
By defining a diagonal matrix as follows 

/ . . . \ 



D := 



6 
6 










• • • U J 



we can write J2i o^Ci — X^DX. This leads us to the following optimization problem 



maximize X'^DX 
s.t X^GX < 1. 



Now, let Y — y/DX. Then, (3-21) can be rewritten as 



maximize Y^Y 

s.t Y'fD-^GD-^Y < 1. 



Suppose D-2GD-2 = G. Then, we have 



maximize Y'^Y 

s.t yt(5y < 1. 



(3-19) 



(3-20) 



(3-21) 



(3-22) 



(3-23) 



The determinant of G—\I = determines the feasible region provided that, A coincides with A'. 
Now, in order to show the ability of our method, we calculate the optimal failure probabilities 
corresponding to unambiguous discrimination of two and three linearly independent states 
with arbitrary prior probabilities. In the simple case of two pure states iV'i) and \ip2) with 
arbitrary prior probabilities 771 and 7)2, G is given by 



G^ 



( ^ 







{ 1 



ai2 



V 







"-12 









1 



/ 



(6) 



ai2 



(3-24) 
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+ - .... 



A2-(^ + ^)A+-V(l-|ai2l) = 0. (3-25) 



Thus, 

66A' - A(6 + 6) + (1 - I«i2l) - 0. (3-26) 
In the feasible region we have pi = A^i and p2 = A^2, then the Eq. (3-26) is equivalent to 

PiP2 - {P1+P2) + (1 - |a?2|) = 0. (3-27) 

This equation determines the feasible region (see Figure. 1). To calculate the minimum 

probability of inconclusive result we put equal the gradient of line rjipi + 772^2 to the gradient 
of equation (3-27), then we will have 

P2 — 1 — Arji and pi — 1 — A772 (3-28) 

By substituting the equation (3-28) into (3-27), we obtain 

(1 + Ar^i) (1 + A7]2) - (2 + A) + (1 - 0^2) = 0, (3-29) 

which implies that A = ±-|^=. Substituting A into equation (3-28), gives the following 
solutions: 

Pi = l±J — |ai2| and p2 = 1 ± W — |ai2|. (3-30) 

V ^1 V 

Since Pi < 1, thus we conclude that 

Pi = 1 - J — |ai2| and P2 = 1 - W — |ai2|. 

V ^1 V 

From the positivity of pi and p2, we have: 

\ai2\<,^< ^ 



In this case, the minimum probability of inconclusive result is: 



(5 = 1- (r/ipi + 772^2) = 2^^1772 1 ai2 1- (3-31) 
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If < |ai2|, then p2 = and if < ^J^, then pi = 0. Then, by using (3-27), we obtain 



Q 


= r]2 + r]i\0'i2? 


■' y VI 


ai2 




' Q 


= 2^771772 |ai2| 


if \ai2\ < y 


fni < ^ 

/ ?7i — |ai2| 


(3-32) 


Q 

V 


= ^1 + ^2|ai2p 


if r^< 


V 





Here, we discuss geometrical interpretation of optimal unambiguous discrimination of two pure 
states on Bloch sphere. One can show that the minimum inconclusive result for unambiguous 
discrimination of two pure states is equivalent to distance between sphere center and the line 
connecting pi to p2 (see Figure 2). Density matrix for a pure qubit state is defined in the 
Bloch form as follows 

Pi = 1{I2 + ni.a) 

P2 = |(/2 +^2.(7). 

For unambiguous discrimination we will have 



(3-33) 



rr(nip2)^o, rr(n2Pi)^o 

Tr(nipi) ^ 1, Tr(n2P2) ^ 1, 
where Hi and 112 are the POVM elements in the pure Bloch form 



(3-34) 



(3-35) 



ni = i(/2 + n'i.(j) 
Tl2^\{l2 + n'^.o). 

This is clear from (3-34) that Tr{Ilip2) is minimal if 71^.712 < 0. Then the optimal case is 
attained for n'^ — —712 and Tr{Jl2Pi) is minimal if n'^ — —rii then the POVM elements are 
given by 

III = lih - n2.a) =1 ijJi){^i I, 

n2 = 1(^2 - ni.a) =1 V'2)('02 I ■ 
For r]\ — r]2 — 1/2, using (3-32) optimum p's coefficients corresponds to 0i = 02 and p's are 
given by 

Pi =P2 = 1 - |ai2| = 1 - cos^, (3-37) 



(3-36) 
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where 29 = ni.n2 and then minimum inconclusive result is equal Q = cos 9. 

Regarding to figure 2 we see that R — cos 9 — Q. That is the minimum inconclusive result 
for unambiguous discrimination of two pure states is equivalent to distance between sphere 
center and the line which connecting pi to p2- 

Now, we give analytical solution for three linearly independent normalized state vectors j'^i) 
, \ilJ2) and ip3 in the three-dimensional complex vector space with arbitrary prior probabilities 
Tji, 7)2 and 773. If we consider — then G is given by 



G 







V 









1 









1 



1 ai2 ai3 



'■12 



1 023 

1 



23 



The characteristic equation is given by 

1 












1 









1 



/ 



«1 



ai2 



6 



1 



''12 



^'23 



13 



;^3 _ _^ j_ ^ — + {— \ ^— 

C2 ^3 C1C2 ^1^3 C2C3 C1C2 C2C3 C1C3 



+ 



1 



/ 2 I 2 

\(l\2 "T 



«13J - 



2012023^13 



1 



Since for feasible region we have pi = A^i, we have 



^166 



ai3 ^ 

_i_ 

S3 / 

(3-38) 



(3-39) 



P1P2P3 - (P2P3 + P1P3 + P1P2) + (1 - 023)^1 



-F(1 - 0^3)^2 + (1 - 012)^3 + a?2 + ^23 + «13 - 2012022013 "1=0 



(3-40) 



If the gradient of the plane 771P1 -|- 772P2 + ?73P3 be equal to the gradient of Eq. (3-40), we will 
have 

f 

P2P3 - {P2 + P3) + 1 - ai3 = Arji 
P1P3 - (Pi + P3) + 1 - 0^3 = Arj2 (3-41) 
P1P2 - {Pi + P2) + 1 - a?2 = ^^3- 
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By solving the Eqs. (3-40) and (3-41), one can obtain 

_ g — «13«23 _ "i2(ai3V^+a23V^) 



/VlV2 



A = 

»73 ' 



A 



ai3a23 I Qii2(ai3v/m+a23V^) 



^ 013 023 _ "l2(ai3Vm-a23v^) 

~ v"7l»?2 y/VlV2V3 

_ _ai3223 I Qi2(ai3V^-a23v/^) 



(3-42) 



but, only A = and A = + "^^^""/^""^^^ give the acceptable values for pi. If the 

point of contact lies in the first octant and < < 1, it gives the optimal solution. If not, 
then an optimal contact point occurs on one of the coordinate planes or even at one of the 
vertices. 

Example 1: We assume that all of the prior probabilities are the same and equal to 1/3. 
If all of the overlaps are the same, i.e., ai2 = Ois = 013 = s where s a is real and positive 
number, then by using equations (3-35) and (3-36) we will obtain 

pis^ + p3S^ + A/3(l - P2) - + = 
Pis2 + p^s^ + A/3(l - ps) - 2s2 + 2s3 = (3-43) 
P2s' + Pss^ + A/3(l - pi) - 2s^ + 2s3 = 

In this case, A = gives the optimal values pi = P2 = P3 — ^ — s and Q = s. 

We can generalize the above example by considering optimal distinguishability for WBE. 
That is, we consider equiangular tight frame WBE sequences (for example special Grassmanian 
frames) [29]. Let {ipi}fLi be an independent frame sequence such that 



{A I '4'j) = s for j 



{ipi I ipj) = 1 for z = j 



(3-44) 



Thus, for optimal distinguishing between N independent vectors which are prepared with equal 
prior probabilities, we can prove (similar to example 1) that, optimal p^'s are given by 



Pi = ... =Piv = 1 - s. 
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One can also prove that, optimal distinguishability corresponds to equal measurement proba- 
bilities p that are equal to the minimum eigenvalue of frame operator. This can be proved by 
defining the frame operator as 

N 

S^Y.\^i){^il (3-45) 

i=l 

such that 

N 

Ski = J2ik\^i){A\l)- (3-46) 

i=l 

Then, we have S = AA^ where A^i = {k \ ipi). On the other hand we define Gramm matrix as 
follows 

G = AU. (3-47) 

One can show that S and G possess equal eigenvalues, thus we find eigenvalues of Gramm 
matrix instead of those of the frame operator. By using (3-45) and (3-48), one can easily see 
that for equiangular tight frame, the Gram matrix can be written as 

G' = / + s(C-l), (3-48) 

where C is a matrix such that all of its matrix entries are equal to one, thus we have — NC. 
Therefore, eigenvalues of G are 

l + s(iV-l), 1-s, (3-49) 

with s < 1. Then, the minimum eigenvalue of the frame operator is equal to 1 — s and thus p 
is equal to 1 — s. 

Example 2: Consider ai2 = 013 = Si and 023 = S2, where both si and S2 are real and 
positive. Then by using equations (3-40) and (3-41) we will have 

P\sl + P3S? + A/3(l - P2) - - s| + 2sls2 = 

Pisl + P2sl + A/3(l - Pa) -si- si + 2sls2 = (3-50) 
P2sl + p^sl + A/3(l - pi) - 2sl + 2sls2 = 0. 
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By solving the above equations we obtain: 

Pi = 1 - 2si 



A = 



^g = 2/3(2si-S2), 

P2 = P3 = S2 - Sl + 1 



= 1/3(^ + 2^2). (3-51) 

P2 = P3 = 1 - S2 



One of the A = 3s^ — 6siS2 and A = which gives smaller inconclusive answer, gives the 
optimal value provided that lies in the feasible region. If is not, then a contact point that 
it is optimal occurs on one of the coordinate plane or even at one of the vertices. Here we 
give two numerical examples such that, in one of them the contact point in equation (3-51) 
is optimal solution, whereas in another one it is not. First, consider the case in which the 
ensemble consists of three linearly independent states with equal prior probabilities 1/3, where 

|^i) = [l,0,Or |^2) = i[l,2,2f |V^3) = |[l,2,-2f. (3-52) 
and si = ^,S2 = |. Then, the optimal solution is given by 

Pi = ^, P2 = Pa = ^ and Q = ^- (3-53) 

In figure 3 we try to show the feasible region of this case which is a convex region. 
As an another example, we consider the case in which the ensemble consists of three state 
vectors with equal probabilities 1/3, where 

l^i) = 73[l'l'ir 1^2) = ^[1,1, of 1^3) = ^[0,1, If. (3-54) 

In this example the contact point occurs on one of the coordinate plane and optimal solution 
is given by: 

pi = 0, p2 = P3 = ^ and Q = ^- (3-55) 
Figure 4 shows the feasible region of this example. 
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3.1 Equal-probability measurement 

A simple measurement that has been employed for unambiguous state discrimination is the 
measurement with pi — p, for all i — 1,2,...,N. This measurement results in equal proba- 
bility of correctly detecting each of the states and is called Equal-probabihty measurement 
(EPM). Using the feasible region, we are able to calculate the prior probabilities, so that EPM 
is optimal. Using equations (3-40) and (3-41), the prior probabilities in the optimal EPM 
measurement for unambitious discrimination of three states are given by 

3 — 2af2 — dis — cii2Cti3 ~ 2a23 + 2023012 + 012 

^ 1 ~ Q12 ~ ^13 "I" Q12Q13 

3 — 2af2 — ^13 — 012^13 ~ 2a23 + 2023^12 + ^12 

773 = 1 - 771 - 772. (3-56) 

In this case the value of p is calculated from equation (3-40). If we consider more than three 
states, the functionality of the 77^ in terms of {aij} are too complicated to be written down, 
so it is not included here. However, if we consider geometrically uniform states, the problem 
will be easy. Let S — — Ui\'^),Ui e ^} be a set of geometrically uniform (GU) states 
generated by a finite group Q of unitary matrices, where jV') is an arbitrary state. Now, let $ 
be the matrix with columns \'4^i). Then, the measurement which minimizes the probability of 
an inconclusive result could be reduced to an equal-probability measurement [27] and consists 
of the measurement operators 

Tli^Pi\i^i){'^i I, (3-57) 

where, Pi is the inverse of the maximum eigenvalue of frame operator [27] and | ipi) = Ui \ 
i^), UieQ with 

I ^) = ((|)(|)*)-i I (3-58) 

In this case, using the feasible region it is easy to show that for optimal EPM measurement 
all of prior probabilities are equal. In general, similar to the example 1, one can prove that for 
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optimal EPM measurement all of pi's are equal to the inverse of maximum eigenvalue of frame 
operator. 

Examples: We now consider an example of a set of GU states. Consider the group G of 
4 unitary matrices Ui, where 



1 ^ 



0-100 
10 

-1 



, Us 



^ 1 ^ 



10 
0-10 
-1 



, C/4 = U^Us. (3-59) 



Now, let the set of GU states is given by S = = Ui\jp),l < i < 4} with = 

1/(3^2) [2, 2, 1, 3]^. Then, we obtain 



3V2 



^ 2 2 2 2 ^ 



2- 2 2-2 
11-1-1 

3- 3-3 3 



(3-60) 



It should be noticed that, the reciprocal states \ipi) = Ui\ip) for i = I, ...,4 with 



^) = ($$*)t I ^) ^ 



4V2 



3 
6 



V2/ 



(3-61) 



are also GU states with generating group Q. Therefore, we can provide the elements of POVM 
as Hi — Pi \ \. Using feasible region it easy to show that the optimal contact point is 

given by 



Pi =P2 =P3 =P4 = g, 

then the equal probability measurement operators are given by 



(3-62) 



(3-63) 



Optimal Unambiguous Discrimination of Quantum States 
where these results are in agreement with those of Ref.[27]. 



20 



4 Lewnstein-Sanpera decomposition (LSD) as an opti- 
mal unambiguous discrimination 

The idea of Refs. [30, 31] is based on the method of subtracting projections on product vectors 
from a given state, that is, for a given density matrix p and any set V — of the states 

belonging to the range of p, one can subtract a density matrix p' = J2iPi^i (not necessarily 
normalized) with Pi>0 such that Sp — p — p' > 0, in the sense that Tr{p') < 1. 

In the following we recall some important definitions and theorems about LSD. 
Definition 1. A non-negative parameter p is called maximal with respect to a (not necessarily 
normalized) density matrix p and the projection operator 11 =| | iff p — pH > 0, and for 

every e > 0, the matrix p — (p + e)n is not positive definite. The maximal p thus determines 
the maximal contribution of 11 that can be subtracted from p maintaining the non-negativity 
of the difference. 

Lemma 1. p is maximal with respect to p and 11 = iff: 

1. if 1^) ^ n{p) thenp = 0, 

2. if IV') e n{p) then 

< p = ~ I ~ . (4-64) 

Definition 2. We say that a pair of non-negative (pi,P2) is maximal 'with, respect to p and a 
pair of projection operators Hi = |'0i)('0i|) n2 = |'02)('02| iff p—pi^i—P2^2 > 0, pi is maximal 
with respect to p — P2n2 and to the projector Hi, p2 is maximal with respect to p — pilli and 
to the projector 112, and the sum pi -|-p2 is maximal. 

Lemma 2. A pair (pi,p2) is maximal with respect to p and a pair of projectors (111,112) iff: 
(a) if l^/^i), \ip2) do not belong to 'R{p) then pi = P2 = 0; (b) if \ipi) does not belong to 7^(p), 
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while 1^2) e 7^(p) then pi = 0, p2 = {H-p\i'2y^; (c) if |V^2) G 7^(p) and (V'i|^|V^2) = 
then Pi = « = 1, 2; (d) finally, if l^-i), |i/'2) e and (V'il^lV'2) 7^ then 



Pi = 



'(^2 1 -1^2) -1(^1 1 -1^2)1), 
P P J 



P2 = ;^((V^l|^IV'l)-|(V'l|^l^2)|), 



(4-65) 
(4-66) 



where D = {i^il^li^i) {H-,\i'2) - \{H-,\i'2)\'. 

Lemma 3. Let a hermitian density matrix p has a decomposition of the form p = p' + (1 —p)6p, 
where p' is a part of density operator p which has the structure p' = Z^^^iPillj, with Ilj being 
projection operator onto state and J2^=iPi — 1- Then the set of {pi}, which are maximal 
with respect to the density matrix p and the set of the projection operators {Hj} form a 
manifold which generically has dimension n — 1 and is determined by the following equation : 



1 - X] ^iPi + 12 DijPiPj - J2 DijkPiPjPk + - + (-1)" PiiPi2-PinD 

i i<j ijk 



ll«2-..ln 



(4-67) 



where the set of {-Dii,i2...im} subdeterminants (minors) of matrix D defined by 



D 



( . . 
an ai2 



Q2I Ct22 



^2n 



(4-68) 



with := (V'j I ^ I "^i)- Equation (4-67) determines feasible region via reciprocal states, that 
is it gives the domain of acceptable values of pi. One way to drive the equation (4-67) is using 
semidefinite programming [15, 17]. 

In the rest of this section, we show that optimal unambiguous discrimination for N linearly 
independent states can be reduced to LSD method. Suppose a quantum system is prepared 
in a state secretly drawn from a known set IV'i), \'^n) where each is a pure state in the 
Hilbert space 7i. In order to discriminate iV'i),..., unambiguously, one can construct a 
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most general POVM consisting of + 1 elements Ho, Hi, IIjv such that 

N 

Ui>0 , i = 0,l,...,N, and ^0^ = /, (4-69) 

i=0 

where / denotes the identity matrix in H. Each element Ilj, i = 1, of POVM corresponds 
to an identification of the corresponding state \ipi), while Ho corresponds to the inconclusive 
answer. For the sake of simphcity, we often specify only IIi, ...,njv for a given POVM since 
the left element Ho is uniquely determined by 

N 

Uo^i-J2^i^i- J2pSi)i^i\- (4-70) 

1=1 i 

The goal of LSD is maximizing p^'s such that Y^iLiPi is maximized. With assuming that 
density matrix in LSD method is equal to identity we obtain LSD as follows 

p-j2^i^i-Y.pSi)m- (4-71) 

i 

With comparison this relation and (4-70) it is clear that we can maximize success probability 
by using LSD method or minimize the inconclusive probability is minimized. Then we say 

that LSD is the same as OptUSDM and we use LSD in order to obtain the elements of the 
optimal POVM. 

4.1 Analytical calculation of optimal POVM for unambiguous dis- 
crimination of quantum states via Lewnstein-Sanpera decom- 
position 

In this section an analytical solution for unambiguous discrimination of two states by using 
Lewnstein-Sanpera decomposition is presented. For three linearly independent states, as the 
LSD method leads to a set of coupled equations which are in general difficult to solve, we 
embark on KKT (see Appendix I) method which makes the problem strongly easy. Since this 
condition is necessary and sufficient, then the answer will be exactly optimal for unambiguous 
discrimination. 
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4.1.1 Optimal unambiguous discrimination of two states 

Suppose that, two pure states iV'i) and IV'2) with arbitrary prior probabihties rji and 772 are 
given. In order to obtain optimal POVM set for these two states by using LSD, we use 
the Lemma 2 of LSD for two states and \ip2)- Let corresponding density matrix in 
Hilbert space 7i is identity operator, i.e., p = I, and II'j^ = rii\ipi) {ipi] and — V2\4'2){4'2\- 
Then, a pair (pi,p2) is maximal with respect to p and the pair of operators Il[ and Ilg iff 
P ~ PiVili^i) {i^i} ~ P2?72|'02)('02| > 0. Therefore, from Lemma 2 of LSD, the maximal pair 



(pi,P2) is given by 



Pi 



a22-y ^|ai2| 
aiia22-|ai2|^ 



(4-72) 



P'^ aiia22-|ai2P ' 

where Ojj = (ipilipj). If the condition 





ai2 


2 


On 


2 + 


ai2 


2 





flu 


2 


an 


2 + 


ai2 


2 



< 771 < I- I ' (4-73) 



is hold, then equation (4-72) is optimal solution. 
If Vi ^ V2, then the optimal solution reads as 



Pi = 0, P2 = ^. (4-74) 

^22 

If ?72 < ?7i; then the optimal solution is given by 

P2 = 0, pi^J-. (4-75) 
an 

Examples 4: For an example we consider following case: Alice gives Bob a qubit repaired in 
one of two states 

1^1) =|0) , 1^2) = ^(1 0)+ I 1)). (4-76) 

Since the states lipi) and \ip2) are non-orthogonal, there is no measurement that can distinguish 
them. In order to obtain optimal POVM, we define dual basis as follows 

1^1) = (I 0)- I 1)) , |^2) = V2|1). (4-77) 
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Therefore, pairs pi and p2 are given by 



Pl=P2 = 



2- 
2 



and finally the elements of optimal POVM are obtained as 



ni = ^^(|o) - |i))((o| - (i|), = (2 - V2)|i)(i 



Ho = 1 - Hi - n2. 



4.1.2 Optimal unambiguous discrimination of three lineeirly independent states 

Now, we consider three linearly independent normalized state vectors ipi , 1P2 and ips with 
arbitrary prior probabilities rji, 772 and 773 in the three-dimensional complex vector space. The 
KKT conditions for unambiguous discrimination of three states are given by 



Pl>0,P2>0,P3>0 

< -P2|^2)(^2| -P3|^3)(^3|)^ = ^(/-Pl|^l)(^l| -P2|^2)(^2| " P3 1^3) (^3 1 ) = 

ZiPi = 0, 2;2P2 = 0, Z3P3 = 0, Zi> 0, 2^2 > 0, 2^3 > 

_ Tr{X\iP,){iP,\)^Zi + r)i, TriX\'4^2){H) ^ Z2 + m,TriX\'4^s){'iP3\) ^ z^ + rj^ 771 > 0, 772 > 0, % > 0. 



(4-78) 



After some calculation, one can obtain 



[(1 -P2a22)(l -^3033) -P2P3|a23r][(l -Plflll) -^2012 








[(1 -piaii)(l -^2022) -piP2|ai2|^][(l -^3033) -^2^23 



[(1 -piaii)(l -^3033) -piP3|ai3r][(l -^2^22) -Piai2 








(4-79) 



With the following conditions 




Pi>0, P2>0, P3>0 



(4-80) 
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the optimal answer is attained as 



Pi = 



P2 



P3 



(022033 ~ 


023 


?) 


- -^1^(021633 - 031023) + (021032 - 022031) 


011(022033 

(011033 - 


Ol3 


023p) - 012(021033 - 031023) + 013(621032 - 022031) 
n - yj(Ol2 033 - 032013) + ^^(011632 - 031012) 


011(022033 
(011022 ~ 


O12 


023 
P) 


2) - 012(021033 - 631623) + 613(621632 - 622631) 
+ ^1^(612623 - 622613) - y|^(^11^23 - 621613) 


ail (5,22033 


— ( 


^23 


2) - 612(621633 - 631623) + 613(621632 - 622631) 



If the conditions (4-82) do not satisfied, then we search another optimal solution. If 





023 


2 


622 


2 + 


623 


2 





O22 


2 


622 


2 + 


623 


2 



is hold, the optimal solation is given by 



Pi =0 

~ 1^22 1 1^33 |- 1^23 P 

|a22|-y^|a23l 



P3 



If 





Ol3 


2 


611 


2 + 


613 


2 



|a22||a33|-|a23r 



<Vi< 





611 


2 


|6ii 


2 + 


613 


2 



is hold, the optimal answer is given by 



I"33|-,^^|ai3l 
^ |aii||a33|-|ai3P 

P2=0 



P3 



If the condition 





012 


2 


|oii 


2 + 


612 


2 



|Sii|-^^|ai3| 
|aii||a33|-|ai3|^ ■ 



<Vi< 





On 


2 


\aii\ 


2 + 


612 


2 



is satisfied, the optimal solution reads as 
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If none of the above conditions is hold, in this case two of pi vanish and the optimal POVM 
will be the same as Von Numan measurement. 

Example 5: In this example we consider reciprocal independent states {ipij^^i such that 

{ipi I ipj) = a for iy^j 

(V^, |V^,) = 1 for i^j, (4-88) 

Thus, for optimal distinguishing of independent vectors that are prepared with equal prob- 
abilities, we minimize the inconclusive result IIq given by 

N 

Uo = I-J2p^\^l'^){A\. (4-89) 

i=l 

Now using LS's theorem and some analytical calculations we show that all of pi's are equal, 
i.e., all of the elements of POVM possess equal probabilities, i.e., 

p,=p, = ...=p, = —l—^. (4-90) 

Similar to GU states, one can prove that the optimal distinguishability corresponds to equal 
measurement probabilities Pi = p for all i such that p is equal to the inverse of maximum 
eigenvalue of the corresponding frame operator. In order to prove this fact, we define frame 
operator as 

TV _ _ 

S = J2\^^){^^l (4-91) 
i=l 

such that 

N 

^« = E(^l^i)(^*IO- (4-92) 
1=1 

Then, S — AA^ with Aki — {k \ On the other hand, the Gramm matrix is defined as 
follows 

G = A'^A. (4-93) 

Again, one can easily show that S and G have equal eigenvalues, thus we evaluate eigenvalues 
of Gramm matrix . The Gramm matrix G can be written as 

G^I + a{C- 1), (4-94) 
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such that its eigenvalues are given by 



l + a{N -1) 



(4-95) 



< 



1 — a, 



where a < 1. Therefore, maximum eigenvalue of frame operator is equal tol-|-a(A'" — 1) and 
thus p is given by (1 + a{N — 

4.2 Discrimination of quantum states using approximated linear 
programming 

As solving the problem analytically is so hard, approximated methods are useful for unam- 
biguous discrimination of linearly independent quantum states. The simplex method is the 
easiest way of solving it. The simplex algorithm is a common algorithm used to solve an op- 
timization problem with a polytope feasible region, such as a linear programming problem. It 
is an improvement over the algorithm to test all feasible solution of the convex feasible region 
and then choose the optimal feasible solution. It does this by moving from one vertex to an 
adjacent vertex, such that the objective function is improved. This algorithm still guarantees 
that the optimal point will be discovered. In addition, only in the worst case scenario all 
vertices will be tested. Here, considering the scope of this paper, a complete treatment of the 
simplex algorithm is unnecessary; for a more complete treatment please refer to any LP text 
such as [32, 33]. In the following, we give some examples and compare the result with ana- 
lytical solution. Suppose that, the quantum system is prepared in one of the two pure states 
l^i) and \ip2) with arbitrary prior probabilities rji and 772. To unambiguously discrimination, 
one should calculate the optimal pi's by solving equation (3 — 27). According the equation 
(3 — 27), the three extremal points 



< 




< 




(4-96) 
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together with the origin (pi = P2 = 0) form a polygon which surrounds the feasible region. 
Since the feasible region is not hnear(see Figure 1) while the polygon is hnear, so we use 
approximated simplex method to find optimal answer. For three linearly independent normal- 
ized state vectors 1-01) , 1-02) and with arbitrary prior probabilities rji, 772 and 773 in the 
three-dimensional complex vector space, there is seven extremal points. These points together 
with the origin form a polygon. The optimal solutions, resulting from the LP method, are 
consistent to that ones obtained analytically in sections 3 and 4 with high accuracy. 

Since obtaining analytical solution for more than three states is difficult, the approximated 
LP method seems to be useful for optimal unambiguous discrimination of linearly indepen- 
dent states. This method, not only enhance the speed of calculation but also is done with high 
precision. 

5 Conclusion 

Here in this work by reducing the theory of the semi-definite programming to a linear program- 
ming one with a feasible region of polygon type which can be solved via simplex method, we 
have been able to obtain optimal measurements to unambiguous discrimination of an arbitrary 
number of pure linearly independent quantum states and using the close connection between 
the Lewenstein-Sanpera decomposition and semi-definite programming, we have been able to 
obtain the optimal positive operator valued measure for some of the well known examples via 
Lewenstein-Sanpera decomposition method. Unambiguous discrimination of mixed states is 
under investigation. 



Appendix I: 
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1 Karush-Kuhn- Tucker (KKT) theorem: 

Assuming that functions gi, hi are differentiable and that strong duahty holds, there exists 
vectors ( E R'^, and y e i?"^, such that the gradient of dual Lagrangian L(x*,C*,y*) = 

f{x*) + J2i Qhiix*) + J2iyi9i{x*) over x vanishes at x* : 

hi{x*) — (primal feasible) 
gi{x*) < (primal feasible) 
y* > (dual feasible) 

y*9^{x*) = 

v/(^*) + EC V hi(x*) + Ey; V 9i(x*) = 0. (i-i) 

i i 

Then x* and {Ci,yi) are primal and dual optimal with zero duality gap. In summary, for any 
convex optimization problem with differentiable objective and constraint functions, the points 
which satisfy the KKT conditions are primal and dual optimal, and have zero duality gap. 
Necessary KKT conditions satisfied by any primal and dual optimal pair and for convex prob- 
lems, KKT conditions are also suffcient. If a convex optimization problem with differentiable 
objective and constraint functions satisfies Slaters condition, then the KKT conditions provide 
necessary and suffcient conditions for optimality: Slaters condition implies that the optimal 
duality gap is zero and the dual optimum is attained, so x is optimal if and only if there are 
{QiVD such that, together with x satisfy the KKT conditions. 

1.1 Slaters condition: 

Suppose X* solves 

minimize f{x)gi{x) >bi,i = l,...,m, (1-2) 
and the feasible set is non empty. Then there is a non-negative vector ^ such that for all x 
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L{x, C) = fix) + C{b - g{x)) < fix*) = L{x\ ()■ 



(1-3) 



In addition, if /(.), gfj(.), i = 1, m are continuously differentiable, then 



dfjx*) 
d{xj) 



-c 



d{x) 



0. 



(1-4) 



In the spatial case the vector x is a solution of the linear program 



minimize 



T 
C X 



s.t Ax — bx>0, 



(1-5) 



if and only if there exist vectors ^ G i?'^, and y G R"^ for which the following conditions hold 



A solution {x*, (*, y*) is called strictly complementary, if x* + y* > 0, i. e., if there exists no 
index i e 1, m such that x* = y* — 0. 
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Figure Captions 

Figure- 1: Feasible region for unambiguous discrimination of two linearly independent 
states with ai2 =< ipi \ 1^2 > is showed by shadow region and approximately feasible region is 
showed by a polygon. 

Figure-2: Unambiguous discrimination of two pure states in Bloch sphere. 

Figure-3: Feasible region for unambiguous discrimination of three linearly independent 
states with ai2 = 013 = | and 023 = ^■ 

Figure-4: Feasible region for unambiguous discrimination of three linearly independent 
with ai2 = ai3 = and 023 = | states. 



